Physical Review Letters, in press 



Heating of the Solar Corona by Dissipative Alfven Solitons 

K. Stasiewicz 
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Solar photo-spheric convection drives myriads of dissipative Alfven solitons (hereinafter called 
alfvenons) capable of accelerating electrons and ions to energies of hundreds of keV and producing 
the X-ray corona. Alfvenons are exact solutions of two-fluid equations for a collisionless plasma and 
represent natural accelerators for conversion of the electromagnetic energy flux driven by convective 
flows into kinetic energy of charged particles in space and astrophysical plasmas. Their properties 
have been experimentally verified in the magnetosphere, where they accelerate auroral electrons to 
tens of keV. 
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PACS numbers: 52.35.Sb, 95.30.Qd, 96.60.Hv, 96.60.Pb 

Understanding the mechanisms that heat plasma in 
the solar corona to temperatures of millions of Kelvin 
has been a long standing problem in solar physics. The 
early view was that the convection in the photosphere 
produces sound waves, internal gravity waves, and mag- 
netohydrodynamic waves which propagate upward to the 
solar corona and deposit their energy to the ambient gas. 
Later, it was realized that all but Alfven waves are dis- 
sipated and/or refracted in the chromosphere and in the 
transition region before reaching the corona 0, 0, LJ U • 
The electromagnetic energy driven by convective flows 
can be transported along the magnetic field as Poynting 
flux of Alfven waves. However, Alfven waves are disin- 
clined to dissipate in collisionless plasmas, and the main 
problem was to explain how this energy flux is deposited 
locally to heat particles in the solar corona |g. Dur- 
ing the last 50 years there have been many attempts to 
solve this outstanding problem in astrophysics, and there 
are more than 20 different models and mechanisms for 
coronal heating proposed in the literature; see reviews 



In this Letter I show that two-fluid equations for a col- 
lisionless plasma have nonlinear solutions in the form of 
dissipative Alfven solitons (alfvenons), which represent 
filamentary space charge structures with strong perpen- 
dicular electric fields and large parallel potential drops. 
Alfvenons will form spontaneously when a magnetohy- 
drodynamic perturbation propagates upward in the solar 
corona where the Alfven speed decreases. They represent 
natural plasma accelerators for converting electromag- 
netic energy flux driven by convective flows to kinetic 
energy of charged particles on spatial scales related to 
the ion inertial length A^ = c/uj p i 1 where c is the speed 
of light and uj p i is ion plasma frequency. Alfvenons can 
provide an explanation for various aspects of electromag- 
netic energy dissipation and heating in the solar corona 
and in the planetary magnetosphcrcs. 

Consider the geometry of magnetic fields in the so- 




photosphere 



FIG. 1: Geometry and coordinate system for a wave struc- 
ture (alfvenon) propagating along x at angle a to the magnetic 
field (blue lines), and driven by convective flows in the pho- 
tosphere. The electric potential structure (red lines) would 
accelerate electrons toward the chromosphere and ions out of 
the chromosphere. The angle a — sin -1 (Vh/Va) is small. 



lar corona depicted in Fig. Q] It shows magnetic field 
lines (blue) inclined at angle a from the wave prop- 
agation direction x. There is a transverse flow with 
speed Vh in the z direction, associated with the elec- 
tric field E = — V/i x B. This geometry implies that 
the photospheric convection drives electromagnetic en- 
ergy flux S = [1q E x B with an upward component S x — 
/Iq 1 Vi 1 B 2 sinacosa. Dragging of magnetic field lines by 
variable convection of a conductive plasma induces mag- 
netic stresses propagating along B with Alfven speed 
Va = B(fiopi)^ 1 ^ 2 , which implies that sin a = (Vhi/Va- 
The upward energy flux carried by Alfven waves driven 
by convection (Vh) is then 
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S x = { P ,/^) ll2 {V h ) 2 Bcosa. 



(1) 



This form of the energy flux is consistent with an em- 
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FIG. 2: Occurrence region of nonlinear waves in two-fluid the- 
ory for plasma /3 = 1CP 2 , with marked position of alfvenons. 
Color coded is log(fcAi), the exponential growth rate of non- 
linear waves. White regions are occupied by linear waves. 



pirical relation between the total energy deposition and 
photospheric magnetic flux in active regions, S x oc B K 
with k rs 1, while k rs 2 should be observed if the en- 
ergy dissipation was related to reconnection of tangled 
magnetic field lines 0|. It is estimated that an average 
upward energy flux of 10 7 erg cm~ 2 s" 1 is sufficient to 
account for energy deposition in the solar corona |lj . For 
the magnetic field B w 100 G and photospheric den- 
sity of pi = 10~ 7 g cm" 3 , a convection speed as small 
as Vh = 0.5 km/s provides Poynting flux exceeding 10 
erg cm -2 s _1 (10 kW m~ 2 ). It is therefore justified to 
assume that the electromagnetic energy flux of Alfven 
modes driven by the photospheric convection (Vh > 0.5 
km/s) is sufficient to power all processes in the chromo- 
sphere and the solar corona. It will be shown below that 
the dissipation mechanism for this Poynting flux can be 
provided by nonlinear Alfven structures similar to those 
observed in association with auroral acceleration. 

The height distributions of the background magnetic 
field B(x) and of ion density pi (x) in the solar atmosphere 
is such that the Alfven speed decreases with x from a 
maximum value inside the chromosphere [8j . It is known 
that Alfven waves can become evanescent or nonlinear 
if they propagate in a region of decreasing Alfven speed 
|lC|. One expects then formation of nonlinear Alfven 
waves (alfvenons) at some altitude in the corona. From 
the general dispersion equation for a two-fluid plasma 
model one can determine the locations of nonlinear waves 
in the phase space (M, a), where M = V x o/Va = w/kVA 
is the alfvenic Mach number. Figure [21 shows portraits 
of the exponential spatial growth rate k\ of nonlinear 
waves 10] . White areas in Fig. [21 correspond to the 
well known linear waves: slow magnetosonic (sM), sound 
(5), Alfven (A), kinetic Alfven waves (kA), fast magne- 
tosonic (fM), and whistler modes (W). The position 
of the sound line (S) is determined by the sound speed 
V s = (-//3/2) 1 / 2 Va, where 7 is the polytropic pressure 



exponent and (3 denotes the ratio of plasma/magnetic 
field pressures. A linear Alfven wave from region A in 
Fig.[2]can tunnel to the nonlinear alfvenon region (below 
the line M = cu/kVA{x) = cos a), if Va(x) is decreasing 
along the wave path. 

In a low beta plasma, details of the pressure model 
|llj | are not relevant and the electron inertia is not im- 
portant for modes propagating quasi-parallel. The gov- 
erning equations for a magnetohydrodynamic structure 
propagating along x can be written in a stationary wave 
frame as Ea IH 
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e x = -M a by tana ^T^TJ - >( 5 ) 

with M a = Mp/cosa, (3 — 2u p /Bq, 6 20 = sina+g(x), 
where g(x) represents a weak gradient of the background 
magnetic field, which can be used to study the transition 
between the linear and nonlinear regimes. Subscript '0' 
denotes background quantities at the starting point xq. 
The magnetic field is normalized with Bq, the electric 
field e x = E x /VaBq, and n — N/Nq is the ion num- 
ber density. Equations ©-© represent a complete set 
of fully nonlinear Hall-MHD equations for field variables 



b y ,b z 



Note that b 7 



cos a, e v — Mq sin a, e z 







are constant, g(xo) = 0, and the variations of the back- 
ground density due to gravitation are neglected. 

The above equations describe linear (sinusoidal) as 
well as nonlinear (cnoidal) waves, including solitons, in 
branches: slow/fast magnetosonic, Alfven, acoustic, and 
kinetic Alfven waves. They are easily integrated with 
an initial perturbation 5b z , implying the boundary val- 
ues: by(x ) = 0, b z {x ) = Sb z , n(x ) = 1 + (M, 2 - 
r y{3/2)~ 1 b z oSb z . Figure shows exact solutions for the 
electric field E x , and the electric potential $ = — J E x dx. 
The magnetic signatures and field-aligned current asso- 
ciated with this structure are shown in Fig. 0J The struc- 
ture has an oppositely directed electric field (V • E < 0) 
with a negative potential in the center, i.e. it is an ion 
hole. The total magnetic field has a small depression 
correlated with the density depression. The computa- 
tions were made in dimensionless variables and converted 
to physical units using typical values measured in the 



B 



10 G and JV k 10 cm , which imply 



Va ~ 2000 km/s, and Xi ~ 20 m. The computed elec- 
tric field can be schematically presented as equipotential 
(red) lines shown in Fig. ^ It is seen that the magnetic 
field lines external to the structure have potential equal 
to zero. The chromospheric potential is also zero because 
it is along the convection streamlines. However, the cen- 
ter magnetic field line is at a large negative potential 
(-80 kV in Fig. 3), implying the existence of a parallel 
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FIG. 3: Exact solutions of equations 12151 for an alfvenon 
propagating at angle a = 0.5°. Medium parameters: Bo = 10 
G, iV = 10 8 cm -3 , (3 = 10" 2 , 7 = 1.6, and A, ; = 20 m. 
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FIG. 4: Field-aligned current density ju (arbitrary units, 
jmin = ^2 x 10~ 8 A m~ 2 ), and transverse magnetic field 
perturbations [b y ,b z , normalized with _Bo) for the alfvenon in 
Fig- El 



potential drop equal to the perpendicular potential drop. 
Conservation of total energy of a particle (kinetic + elec- 
tric potential) implies that an electron leaving the struc- 
ture along the magnetic held lines to the chromosphere 
will acquire kinetic energy equivalent to e$(x), produc- 
ing X-rays when interacting with the ambient gas. Sim- 
ilarly, an ion from the chromosphere entering the struc- 
ture would be accelerated upward by the same potential 
$(x). Mathematically, this nonlinear structure is formed 
as a result of a balance between nonlinear growth and 
dispersion related to terms with ion inertia Xi in the gov- 
erning equations ©-©, and therefore it could be re- 
garded as a soliton. However, connection with the con- 
ducting chromosphere introduces a dissipative element 
to this soliton. Its energy would dissipate by accelera- 
tion of electrons toward the chromosphere and ions out 
of the chromosphere. Alfvenons form filamentary struc- 



tures with extension along the magnetic field much larger 
than the perpendicular size. In the example of Fig. 3, 
the length of the alfvenon is L ~ 3 x 10 5 Aj ~ 6000 km, 
while the width L± sw Lsina « 50 km. The length of 
the alfvenon would increase (decrease) with decreasing 
(increasing) angle a = sin (Vh/Va)- 

Let us estimate the electric field available for accel- 
eration of charged particles in alfvenons. From Eq. (J5J 
we find that E x w — MVaB v tana. A typical magnetic 
polarization pattern for alfvenons, which can be found 
from numerical solutions, is rotation of the transverse 
magnetic field (B y ,B z ) around the guiding B x field such 
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B z q = Bq sin a. Furthermore, alfvenons 



are formed near M « cos a (see Fig- EJ), which implies 



\E X 



VaB sin a. 



(6) 



As is seen, the electric field in an alfvenon depends on 
the propagation angle and can vary from zero (parallel 
propagation) to the maximum value of Ea = BqVa for 
perpendicular propagation. 

The electric potential structure obtained here in a self- 
consistent way (though in a simplified geometry) is not 
a speculative result awaiting observational verification. 
In fact, the presence of such electric potential structures 
above the aurora was first inferred from satellite measure- 
ments 30 years ago [l3| . Properties of these U-shaped au- 
roral acceleration structures have been thoroughly inves- 
tigated in numerous publications listed in a recent review 
of auroral processes [l4| , but until now there was no the- 
oretical model for these structures. Thus, the alfvenons 
described in this Letter provide also an explanation for 
the U-shaped auroral potential structures. Particle mea- 
surements inside these acceleration structures show up- 
ward directed ion beams associated with electrons accel- 
erated downward by the electric potentials 1-20 kV, and 
forming so-called inverted-V structures in energy versus 
time flux-spectrograms. An obvious difference between 
the solar corona and the magnetosphcre is that in the 
first case the driving convection is applied to the bot- 
tom end of the flux tubes in Fig. ^ (in the photosphere) , 
while in the second case to the top end (in the distant 
magnetosphere) . The magnetospheric alfvenons are cre- 
ated below an altitude of 1 Re, i.e. in region where the 
Alfven speed starts to decrease toward the ionosphere, as 
predicted by this theory. Similar acceleration structures 
have been recently measured on Mars |15| , indicating the 
universal applicability of the alfvenon mechanism. One 
should also mention that the present theory of nonlin- 
ear waves, when extended with anisotropic and polybaric 
pressure equations jllj . gives numerical results in quan- 
titative agreement with observations of large amplitude 
(SB /B ~ 200%) trains of magnetosonic solitons in a hot 
(j3 ~ 10) magnetosheath plasma [lflllg- Plasmas in the 
solar atmosphere and the magnetosphere are quite simi- 
lar, with the same range of plasma beta (10~ 3 — 10 1 ), and 
a similar altitude transition from a partly ionized, colli- 
sion dominated to a fully ionized, collisionless plasma. 

The numerical results obtained in this paper indicate 



that depending on the propagation angle and parameters 
Va, B, (3 in different regions of the corona one could pos- 
sibly create acceleration structures with electric potential 
drops that would account for most of the observed X-ray 
and radio emissions due to accelerated electrons. Typi- 
cal integrated voltages across the alfvenon structures in 
the solar corona determined numerically are hundreds 
of kilovolt. Because the photospheric convection is a 
variable and permanently occurring process, myriads of 
alfvenons will be continuously recreated at different al- 
titudes over the whole Sun, producing what is observed 
as the X-ray corona. Alfvenons could form bundles of 
threads or sheets and build up larger current structures, 
which should be seen in X-ray emissions when acceler- 
ated electrons thermalize in denser regions. Actually, 
high-resolution images from TRACE (Transition Region 
and Coronal Explorer) are suggestive of the existence of 
such threads, sheets and arcades of energized particles in 
the solar corona. 

The present theory provides a natural link to explana- 
tions for the acceleration of the solar wind and the cre- 
ation of solar flares. The solar wind originates from the 
coronal regions cool in X-rays, which have open magnetic 
field lines. Solar wind ions have energy concentrated in 
the drift velocity with a small thermal spread (Vd 3> Vt), 
which is consistent with acceleration by electric potential 
difference, and not by any stochastic or turbulent process 
that would produce Vd ~ Vt. As is seen in Fig. Q] an el- 
ementary alfvenon event ejects ions outward, providing 
the first kick-off for the solar wind. Typical speeds of the 
solar wind protons 400-800 km/s could be produced by 
an electric potential of 1-4 kV. A possible explanation for 
such small values of the acceleration voltage in coronal 
holes is that plasma on open field lines cannot support 
significant parallel electric fields. 

On closed coronal loops, periodic energy releases by 
alfvenons driven by photospheric convection would heat 
plasma that will remain confined in magnetic traps. 
This would allow large pressure gradients and pitch-angle 
anisotropy to build up and support increasingly larger 
parallel electric fields that would accelerate ions upward. 
When the magnetic field tension at the top of the loop 
will not be able to balance the Reynolds stresses associ- 
ated with flows of accelerated ions, the loop would break 



up, leading to plasma jets, coronal mass ejections, and a 
topological reconfiguration of the magnetic field. A de- 
tailed description of the mechanism for flares implied by 
the present theory is beyond the scope of this Letter and 
will be addressed elsewhere. 

Application of the present model to heating of the solar 
corona can be summarized as follows: 

(a) Dragging of magnetic field lines by non-stationary 
photospheric convection Vh induces magnetohydrody- 
namic perturbations propagating upward and carrying 
energy flux given by Eq. (|TJ. 

(b) An upward propagating magnetic perturbation be- 
comes nonlinear in regions of decreasing Va (x) and forms 
an alfvenon, i.e. the U-shaped electric potential structure 
shown in Fig. Q] The size of the structures scales with 
Aj and the propagation angle. For a propagation angle 
a w 0.5° and Aj ~ 20 m the size is: Ln ~ 6000 km, 
L± ~ 50 km, and the generated voltage is <& ~ 100 kV. 

(c) The electric potential structures created in the solar 
corona dissipate electromagnetic energy through direct 
acceleration of electrons toward the chromosphere and 
ions out of the chromosphere. 

(d) Ions accelerated outward by 1-4 kV potential drops 
in coronal holes would initiate outflow of the solar wind. 

(e) Electrons accelerated toward the chromosphere by 
a potential difference up to hundreds kV (on closed loops) 
produce X-ray corona, radio emissions, and evapora- 
tion/heating of the chromosphere. 

(f) The processes described in (a)-(e) would repeat as 
long as there is variable photospheric convection pumping 
the energy flux (1) into regions of a decreasing Alfven 
speed. 

The alfvenons introduced in this Letter appear to be 
effective and spectacular converters of electromagnetic 
energy flux into kinetic energy of particles. They have 
been measured by numerous spacecraft in the terrestrial 
magnetosphere, where they accelerate auroral electrons 
toward the ionosphere to tens of keV, and are detected 
also in the martian environment. They appear to be of 
universal importance for astrophysical plasmas and must 
occur also in the solar atmosphere, where they can ac- 
count for the acceleration and heating of plasma in the 
solar corona. 
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